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We derive finite temperature time-dependent effective actions for the phase of the pairing field, 
which are appropriate for a 2D electron system with both non-retarded d- and s-wave attraction. 
As for s-wave pairing the d-wave effective action contains terms with Landau damping, but their 
structure appears to be different from the s-wave case due to the fact that the Landau damping is 
determined by the quasiparticle group velocity v 9 , which for d-wave pairing does not have the same 
direction as the non-interacting Fermi velocity vf. We show that for d-wave pairing the Landau 
term has a linear low temperature dependence and in contrast to the s-wave case are important for 
all finite temperatures. A possible experimental observation of the phase excitations is discussed. 
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The microscopic derivation of the effective time-dependent Ginzburg-Landau (GL) theory continues to attract 
attention since an early paper by Abrahams and Tsuneto 0]. Whereas the static GL potential was derived j| from 
the microscopic BCS theory soon after its introduction, the time-dependent GL theory is still a subject of interest (see 
||, Q for a review on the problem's history). One of the reasons for this is the presence of Landau damping terms in 
the effective action. For s-wave superconductivity these terms are singular at the origin of energy-momentum space, 
' and consequently they cannot be expanded as a Taylor series about the origin. In other words, these terms do not 
have a well-defined expansion in terms of space and time derivatives of the ordering field and therefore they cannot 
7— I ' be represented as a part of a local Lagrangian. We recall that at T — and for the static (time-independent) case 
the Landau damping vanishes, so that either at T — one still has a local well-defined time-dependent GL theory or 
for T =/= the familiar static GL theory exists. It is known, however, that for s-wave superconductivity even though 
the Landau terms do exist, they appear to be small compared to the main terms of the effective action in the large 
temperature region < T < 0.6T C Q, where T c is the superconducting transition temperature. This is evidently 
related to the fact that only thermally excited quasiparticles contribute to the Landau damping. The number of such 
quasiparticles at low temperatures appears to be a small fraction of the total charge carriers number in the s-wave 
superconductor due to the nonzero superconducting gap A s which opens over all directions on the Fermi surface. 
For a d-wave superconductor there are four Dirac points (nodes) where the superconducting gap A^(k) becomes 
I ■ zero on the Fermi surface. The presence of the nodes increases significantly the number of the thermally excited 
quasiparticles at given temperature T comparing to the s-wave case. Therefore one can expect that the Landau 
damping should be stronger for superconductors with a d-wave gap which is commonly accepted to be the case of high- 
temperature superconductors (HTSC) ||. Moreover, it is believed that at temperatures T -C T c , these quasiparticles 
are reasonably well described by the Landau quasiparticles, even though such an approach fails in these materials at 
higher energies Q. This is the reason why one can hope that a generalization of the BCS-like approach |j| for the 
2D d-wave superconductivity may be relevant to the description of the low-temperature time-dependent GL theory 
in HTSC. 

In this work we derive such a theory from a microscopical model with d-wave pairing extending the approach of 
developed for s-wave superconductivity. As known from the physical origin of the Landau damping is a scattering 
of the thermally excited quasiparticles ( "normal" fluid) with group velocity v 5 from the excitations of phase (or 0-) 
quantums. Such conversion occurs only if the Cerenkov condition, £1 = v g K for the energy fl and momentum K of 
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the ^-excitation is satisfied. This phenomenon in superconductivity is also called Landau damping since its equivalent 
in the plasma theory was originally obtained by Landau (see e.g. Q). 

To emphasize the difference between the Landau damping for s- and d-wave pairing, we also derive for the compar- 
ison the corresponding terms for a 2D s-wave superconductor. In addition we compare 2D expressions obtained here 
with the 3D s-wave case studied in M. The collective phase oscillations in charged d-wave superconductor for clean 
and dirty cases were recently studied in |^] . Due to the complexity of the corresponding equations they were solved 
numerically neglecting damping of the phase excitations. Thus our fully analytical treatment can be very useful for 
further studies of the phase excitations. We also mention a recent paper || where the effective action for the phase 
mode in the d-wave superconductor was obtained using the cumulant expansion. The Landau terms were neglected 
in j^], but the effect of Coulomb interaction was taken into account. 

Our main results can be summarized as follows. 

1. We find that the main physical difference between the s- and d-wave cases is related to the fact that for d-wave 
superconductivity the direction of the quasiparticle group velocity v g (k) = dE(k)/dk {E(k) is the quasiparticle 
dispersion law) does not coincide with the Fermi velocity vf || and a gap velocity va = 9Ad(k)/<9k also enters 
into the Cerenkov condition along with Vf . 

2. We show that the intensity of the Landau damping has a linear temperature dependence at low T with a 
coefficient expressed in terms of the anisotropy ao = Vf/va of the Dirac spectrum -E(k) = \J v'pk'f + v^k^. Here 
k\ (k 2 ) are the projections of the quasiparticle momentum on the directions perpendicular (parallel) to the Fermi 
surface. The parameters vp, va and ap, proved to be very convenient both in the theory, for example, of the transport 
phenomena Q , ultrasonic attenuation jlO| in d-wave superconductors and for the analysis of various experiments . 

3. We find that the Landau damping is sensitive to the direction of the phason momentum. In particular, for a 
given node the Landau damping is possible only if the components of the phason momentum K = (K\, K 2 ) (which are 
defined exactly as the components of the quasiparticle momentum above) satisfy the condition \Q\ < yJv F K( + ujifj. 

4. We derive a simple approximate representation for the Landau damping terms which can be useful for further 
studies of the d-wave superconductors. 

5. We derive an approximate expression for the propagator of the Bogolyubov- Anderson mode which includes the 
Landau damping. 

6. Concerning the mathematical formalism used in the paper, we adapt the bilocal Hubbard-Stratonovich field 
method of |]l2| to the d-wave pairing. Additionally we adjust the technique of the derivative expansion for the "phase 
only" action (see Q and Refs. therein) for the model with the tight-binding spectrum. 

The paper is organized as follows: In Section II we present our model and write down the partition function using 
a bilocal Hubbard-Stratonovich field. In Section III, we introduce the modulus-phase variables and represent the 



effective "phase only" action as an infinite series. It appears that the low-energy phase dynamics is contained in the 
first two terms which are evaluated, respectively, in Sections [Fv| and [v] with some details considered in Appendix [a|. 
The eff ectiv e Lagrangians for the d- and s-wave cases without the Landau damping are discussed in Section |Vl| In 
Section VII we derive the damping terms and in detail compare d- and s-wave cases. The approximate forms of the 
effective action and 9- propagator are considered in Section |VII] . Section IX presents our conclusions and comments 
on a possible experimental observation of the phase excitations. 



II. MODEL 



Let us consider the following action 

r0 



r=(x,y), (3=- 



S =-l dT E f d 2 r^(T,r)d T MT,r)+H(r) 
Jo I a J 

where the Hamiltonian H (r) is 

^( T )=E / d 2 ^(T,r)[e(-zV)- M ]Vv(T,r) 

-^E / d2ri J d2r ^l(T,r 2 )'ipU T , r iW( r ii r 2)'ips(T, ri)Vv(i", r 2 ) 



(1) 



(2) 



Here i> a ij-, r) is a fermion field with the spin a =f , 4, a = —cr, t is the imaginary time and V(r\] r 2 ) is an attractive 
potential. For the sake of the simplicity we consider the dispersion law, e(k) = — 2£(cos k x a + cos k y a) for a model on 
a square lattice with the constant a including the nearest-neighbor hopping t only. This, however, is not an essential 



3 



restriction because the final results for the d-wave case will be formulated in terms of the non-interacting Fermi 
velocity vp = 9e(k)/9k|k = k F and the gap velocity va defined in the Introduction. 

The bilocal Hubbard- Stratonovich fields <I>(t, i"i;r 2 ) and &(t, ri; r 2 ) (see e.g. Q) can be utilized to study the 
model ©, |) 



exp 



dr / d 2 n / d 2 r 2 ^ T t (r,r 2 )V'|(T,ri)t/(ri;r 2 )V'i(r,ri)V' T (r,r 2 ) 



= J P$ t (r,r 1 ;r 2 )I?$(T,r 1 ;r 2 )exp - J ' dr J d 2 n J ^ * — |$(r, n; r 2 )| 2 
+ [ dr [ d 2 r x f d 2 r2($ t (T,r 1 ;r 2 )^(T,r 1 )^ T (T,r2)+^(T,ri)^J(T,r2)$(r,r 1 ;r 2 )) 



(3) 



On the right hand side, 1/V^(ri;r 2 ) is understood as numeric division, no matrix inversion being implied. The 
hcrmitian conjugate of <£>(t, ri;r 2 ) includes the transpose in the functional sense, i.e. $^(r, ri; r 2 ) = [$(t, r 2 ; ri)]*. 
Thus in the Nambu variables 



*(t, r) 

the partition function can be written as 

Z = / D* t I?4'I?$ t I?$ exp 



(4) 



1 



V(ri;r 2 ) 



|$(r,r i; r 2 )| 2 



+ *t( T> n)(-a T - r 3 e(-iT 3 V))*(r, r 2 )*(ri - r 2 ) 

+$ + (r, n;r 2 )* + (r, ri)r_tf (r, r 2 ) + *t(r, ri)r + *(r, r 2 )$(r, n; r 2 )] } 



(5) 



where £(— iT3V) = e(— iraV) — /i and T3, t± = (t\ ± ir 2 )/2 are Pauli matrices. 

In general an electron-electron attraction on the nearest-neighbor lattice sites can be considered (see, for example, 
fl3| , |l^|). The momentum representation for this interaction contain in the pairing channel extended s-, d- and even 
p-wave pairing terms: 



V(k - k') = V[cos(k x - k' x )a + cos(k y - k' y )a] 

= — (cos k x a + cosk y a)(cosk' x a + cos k' y a) + — (cos k x a — cos k y a)(cos k' x a — cos k' y a) 
+ V (sin k x sin k' x + sin k y sin k' ) . 



(0) 



Motivating by HTSC we consider here d-wave pairing only, so that for the Fourier transform of the pairing potential 
V(ri — r 2 ) we use 



y(k — k') = Vd(cos k x a — cos fcj,a)(cos k' x a — cos k'a) . 



(7) 



As was mentioned in Introduction, we will compare the main results for the c?-wave case with the simplest 2D 
continuum s-wave pairing model (see e.g. [ fl5[ and the review |l6||) which has a quadratic dispersion law e(k) = k 2 /2m 
and a local attraction V(t\ — r 2 ) = V8(r\ — r 2 ). 



III. THE EFFECTIVE ACTION 

While for the model with the local four-fermion attraction ^, [l5| the modulus-phase variables could be introduced 
exactly, one should apply an additional approximation for treating the present model. 

Let us split the charged fermi-fields ip(r, r) and ip^(r, r) in (||) into the neutral fermi-field x( T : r ) 10 an d charged 
Bose-field exp[i6(r, r)/2] 

Mr, r) = X*(r, r) exp[i0(r, r)/2] , ^(r, r) = xt(r, r) exp[~z0(r, r)/2] . (8) 

It is clear that the terms containing S(ri — r 2 ) in (||) can be treated similarly to the old model and the problem arises 
when one deals with the Hubbard-Stratonovich field. To consider this field we introduce the relative r = ri r 2 
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and center of mass coordinates R = (ri + r 2 )/2. Now we can introduce the modulus-phase representation for the 
Hubbard-Stratonovich field 



$(T,n,r 2 ) = $(r,R,r) = A(r, R, r) exp[i0(r, R, r)] 



(9) 



where A(R, r) is the modulus of the Hubbard-Stratonovich field and 0(R,r) is its phase. 

Assuming that the global phase 0(R, r) varies slowly over distances on the order of a Cooper pair size and thus is 
not sensitive to the inner pair structure described by the relative variable r, we can rewrite (H) as 



$(r, R, r) fa A(t, R, r) exp[i6»(r, R)] 



(10) 



The approximation we made writing Eq. (10) is in fact equivalent to the Born-Oppenheimer approximation [la] 
which allows one to separate the dynamics of the Cooper pair formation described by the relative coordinate r in 
A(r, R, r) from the motion of the superconducting condensate described by the center mass coordinate R in 0(r, R). 
If the condensate motion is slow enough this separation becomes possible because the dynamics of the Cooper pai r 
formation can always follow the motion of the condensate. Using the lattice language one can also say about ([H]) 
that the bond phase is replaced bythe site phase ]l9| . 

Applying the transformations ( J10| ) and (g) to the terms with the Hubbard-Stratonovich field in (||) we obtain (the 
imaginary time r is omitted) 



$t(R;r)^i(ri)V>T(r2) - $ f (R; r) Xi (n) exp 
w A(R; r) exp[-i0(R)]x4.(ri)x T (r 2 ) exp 



(R + r/2) 



R) 



2 

2 2 



Xt(r 2 )exp 
V Q V /3 0(R 



z0(R- r/2) 



(11) 



A(R;r)xi(r 1 )xT(r 2 ) 



where we used the assumption (or hydrodynamical, long wavelength approximation, see [2Q] ) that 0(R) varies slowly, 
0(R) S> £q(V0(R)) 2 . Here £0 is the coherence length which for the BCS theory coincides with an average pair size. 
Then the partition function in the modulus-phase variables is 



Z = J AX>AX>0exp[-/3ft(A,d0)] 



where the effective potential 



/3£)(A, 86) = £ dr J d 2 n J 



2 A (t, R, r) , 
d 2 r 2 T , K ' ' I - TrLnG" 1 
V{ri - r 2 ) 



with 



G _1 = G^ 1 - E , 



(12) 



(13) 



(14) 



£ 1 (Ti,T 2 ;ri,r 2 ) = (n,ri|^ 1 \r 2 ,r 2 ) 

-Id Tl - r 3 ^(-ir 3 Vi) <5(ri - r 2 )<5(ri - r 2 ) + nA(n - r 2 ,R,r) ; 



(15) 



(ri,ri|E|T 2 ,r 2 ) = 
iia 2 V2 ,0 



T3 1 



Or 



ta 2 ^^-coa(-iaV Xl ) + (x -> y) 



-I - 



2 4 
cos(— iaV xl ) + taVa; 1 0sin(— iaV Xl ) + (% — > y) 



<S(ti - r 2 )(5(ri - r 2 ) . 



(16) 



Thus the gauge transformation (g) resulted in the separation of the dependences on A and 0, viz. is present only 
in E. The similar method of the derivative expansion was used before in j|, As pointed out in B the method 
allows to maintain explicitly the Galilean invariance (the Landau terms break it) and the continuity equation, while 
the expansion $(x) = A + $i(x) +«$ 2 (x) used recently in || demands the additional enforcement of the conservation 
laws p2[ . 

Since the low-energy dynamics in the phase in which A 7^ is determined by the long-wavelength fluctuations 
of 0(x), only the lowest order derivatives of the phase such as V0, d T 9 and V 2 need be retained in what follows. 
However, to take into account the tight-binding electron spectrum the operators sin(— iaV) and cos(— iaSI) must be 
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kept. Thus in E we have omitted higher order terms in V0, but in order to keep all relevant terms in the expansion 
of sin(— ietV) the necessary resummation was done One can easily see that for the quadratic dispersion law 

2t — » l/(ma 2 ), cos(— idV) — > 1 and sin(— iaV) — > — iaV, so that Eq. ( |l6| ) reduces to the known expression from 
0, Elf. Thus we arrive at the one- loop effective action 



Q ~ f} fei „ (v, fi, T, A, 30) + Sl%*(v, (j,, T, A) 



where 



oo 1 

n Hn (li,T,A,dd) = TTr V -(£?£)' 



SA/9R=0 



and 



^ F t (M,r,A) 



d 2 i? / d 



A 2 (r) 



- TTr Ln£" 



(17) 



(18) 



(19) 



9A/aR=0 



Deriving the "phase only" action for the s-wave model it was possible to use A(R, r) = const The d-wave case 
is more complicated because one should keep the dependence on the relative coordinate, A(R, r) A(r) which is 
related to the nontrivial pairing. The dependences of the gap A on T, \i and Vd follow from the extremum condition 
for the mean-field (OA/OR = 0) potential t?Q"J t /<3A = which results in the usual BCS gap equation. For the d-wave 
pairing potential (^) one obtains 



Ad(k) = ^^(cos k x a — cos k y a) 



(20) 



where A^ is the gap amplitude. In our case there is no need to solve the gap equation and express Ad in terms of T, 
/j, and V c i since in what follows we will use A^, or more precisely the velocity va, as the input parameters and will be 
interested in the low temperature (T <C Ad) regime. 

Thus assuming that A(R, r) does not depend on R one obtains for the frequency- momentum representation of (|l5| ) 



^ n J + r 3 g(k) -nA(k) 

W 2+£2(k) + A2(k) 



(21) 



where A(k) is given by ( p0| ) and u) n = 7r(2n + 1)T is fermionic (odd) Matsubara frequency. 

The phase dynamics is contained in the kinetic part !l|; m of the effective action which only involves the single 
degree of freedom 9. As discussed in it is enough to restrict ourselves to terms with n = 1, 2 in the infinite series 
in ([l8]) since at T = this would give the right answer for a local time-dependent GL functional which involves the 
derivatives not higher than (V0) 4 and (pt9) 2 . 

IV. THE FIRST ORDER TERM OF THE EFFECTIVE ACTION AND THE NODAL APPROXIMATION 



In this section we calculate the first (n — 1) term of the sum appearing in (18) 
n& = TTr [<?!!] 



= T fdrf d 2 r\ T £ f tr^.k) 

n— — oo ^ ' 

Summing over Matsubara frequencies, one obtains 



i 



d T 9 to 2 



ta z 



H (y x 0) cos k x a H (V y 9) cos k y a 



=T I dr I d 2 r 







T2 n(k) z— + m xx (k) — 



(2tt) 2 



m-i(k)^ 



with m-!(k) = d 2 £(k)/d£; 2 , m^k) = ^(k)/^ 2 and 



»(k) = 1 - |^ tanh < ^ = Ve 2 (k) + A 2 (k) 



(22) 



(23) 



(24) 
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For T <C Arf linearizing the quasiparticle spectrum about the nodes and defining a coordinate system (ki,k 2 ) at 
each node with ki (k 2 ) perpendicular (parallel) to the Fermi surface, we can replace the momentum integration 
in (^3|) by an integral over the k-space area surrounding each node [0]. If we further define a scaled momentum 
P = (piiPz) = (P; v) we can let 



d 2 k 
(2tt) 2 



E 

i=i 



dkidk2 



E 



d 2 p 



(2tt) 2 vfva 



E 



pdp 



2ttvfva 



2 71 



dip 
2^ 



7T" » 



y/TTV F VA , (25) 



where pi = Ui?/ci = e(k) = pcos^, p 2 = «a&2 = Ad(k) = psin<^ and p = \fp\ 



Note that for the particular square lattice model used above those velocities are vf — 2\[2ta and va = A<ja/\/2 
respectively. 

Using ( P5|) one can express Eq. (|2^) in terms of vf and ua 



(i) _ 



T / dr d 



in } - 



d T 9 



2tt 



pdp 

2itvfva J 2ir 



dw p a 2 p cos 2 ip „, 2 

— tanh i^iii- 



2T 



T J dr j <l-r 



in f - 



d T 9 ^ttvfva 
+ 48a 



(V0) 5 



(26) 



where 



rf 2 fc 
(2^) 2 



>(k) 



(27) 



is the density of carriers. We note that due to the slow convergence of the integrand in (j2(j) the final expression 
depends explicitly on the value of the momentum cutoff p m ax which was defined in j(| in such a way that the area 
of the new integration region over 4 Brillouin sub-zones (see Fig. is the same as that of the original Brillouin 
zone. Going from Eq. ( p3[ ) to (26) we essentially replace the averaging over the true Fermi surface of the system by 
the averaging over 4 nodal sub-zones. The validity of this approximation can only be justified if the corresponding 
integrals contain the derivative of the Fermi distribution np(k-) which is highly peaked in the vicinity of the nodes. 
This appears to be the case of the temperature dependent parts of the phase stiffness J(T), compressibility K(T) 
and the Landau damping terms. For the zero te mper ature values J(T = 0) and K{T = 0) the nodal approximation 
is not well justified. However, as we show in Sec. |V4 this approximation can be justified a forteriory for their ratio 
(see Fig. ||) which determines the velocity of the Bogolyubov-Anderson-Goldstone mode. Finally, we stress that after 
approximation is used, it is impossible to recover the s-wave limit by putting va — * 0. 



V. THE SECOND ORDER TERM OF THE EFFECTIVE ACTION 



Let us evaluate the trace of the second term in expansion (h 



^ = |Tr[c;vx;-j 



Substituting ( jlq ) into ( pq ) we obtain that 



(28) 



(29) 



where using {6 ■ ■ ■ 9} we denoted symbolically that the corresponding term of (28) is either diagonal (i.e. its frequency- 
momentum representation contains 9{iVt n , K)f2 2 6(— ifl n , ~ K) or 8(ifl n , K)K^9(— ii7„, — K)) as 



^{^10} = - J2 J j^8(i^,K)\--f)0(-iSl n ,-K)T j^^(^,K; !Wi ,k); 



(30) 



pn^ n {6K 2 e} = | f; j0p 9(iSl n , K)K 2 J(-iQ n , -K) x 



^E 



d 2 k 

WT 2 



(31) 



7r o(i^„,K;iw;,k)(ta) 2 sin(k x - K x /2)asia(k x + K x /2)a+ (x — > y) 
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or mixed 



(3O<»{ffKn n 0} = -T jr j |^ 9(iQ n , K)^9(-iQ n , -K) 



T E 



d 2 k 
(2nf 



(2n) 2 

2,t& zt& 

7ro3(^n,K;iwz,k)— sin(fc a; + K x /2)a + TT 30 (iQ n ,K; iw/,k) — sm(k x - K x /2)a 



(32) 



+ (x -> y) . 

In (|3(]) - ( |32| ) we introduced the following short-hand notations 

7Ty(ifi n , K;iw;,k) = tr[S(icjj + iSl T1 , k + K/2)TiG(i(jji, k - K/2)r 3 -] , = (t = I,t 3 ) . 

More generally, we can rewrite Eqs. ([5l|) and (^) as follows 

n— — 00 



(33) 



(34) 



T 



E 



n— — 00 



6{in n ,K.)- 



ici n , — k) [n«(in n ,K) + ng (ifi»,K)] 



(35) 



where 



00 f d 2 k 

-77 7Too(i^n, K; IW;, k)uF a (k)VFp (k) , 

^ — — 00 

~ r d 2 k 

H% 3 (iQn,K) = 2J / 7To3(^n,K;io;i,k)i; FQ (k) , 



(36) 



«Fa(k) = d£(k)/dk a , we used the approximation u^ Q (k) ~ Up a (k ± K/2) and took into account that 
(l/27r) J d 2 kv a vp = J kdkv 2 8 a p/2. It is convenient to introduce here 



~ f d 2 k 

U 33 (ifl n ,K) = 22 / (2 )2 7T33(^n>K;io;;,k) 
£ — — 00 



(37) 



which would allow to rewrite ( J30| ) in the same fashion as (|3J). 

As one could notice the product of the Fermi velocities enters Eq. ( |34| ) via ( |36| ) . There is nothing surprising in this 
fact since this piece of the effective action is related to the paramagnetic current correlator (j a jp) [§] arid in its turn 
the current operator, j contains the Fermi velocity v^(k). We will return to this point considering the Landau term 
which originates from (|34j ), so that here we note only that the current correlator term along with the diamagnetic 
term ~ (X79) 2 in Eq. fl22|) form together the mean-field phase stiffness. 

The matrix traces ntj and the corresponding expressions for Iljj are calculated in Appendix [A|. Although the 
expressions for them are rather lengthy they have a clear physical interpretation which is also discussed in the 
Appendix. 

VI. THE EFFECTIVE LAGRANGIAN ATT/O WITHOUT THE LANDAU TERMS 



The contribution of the first order term to the effective action is given in (|26|). Concerning the second order term, 
we note that when the Landau terms are neglected it is enough to set f2 n = in side II in (|30|), ( |34"| ) and (|35|). To 
be more precise, the Landau terms arise from the second line of Eqs. (A2) and (A3) which contains "dangerous" 
denominators 1/(E+ — £L ± if2 n ). One can however notice that the second line of Eq. (A2) leads also to the regular 
terms which are proportional to the derivative dnp(E)/dE. These lead to the second term in the square brackets in 
( p8[ ) and the whole expression ( |39| ) shown below. For the s-wave superconductivity |l5| in the temperature region 
< T < 0.6T C these terms are very small compared to the main terms. Although the contribution from these terms 
is still local, their presence breaks the Galilean invariance B . This is the reason why it was more natural for to 
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t rea t them along with "true" Landau terms since they also originate from the same denominators of the second line of 
( A2) as was mentioned above. For the d-wave superconductivity this splitting, however, appears to be rather artificial 
since these terms are not small even for low temperatures due to the presence of the nodal quasiparticles, so here we 
will consider all regular terms. 

For the regular terms from the second order term we obtain a local effective action, involving time and space 
derivatives of 9(t,r): 



= i- 



dt 



d 2 r 



1 f dfL 

2 J 

(d t 6(t,r)) 2 



d 2 K Q 2 
(2ir) 2 y ' ' 4 



d 2 k 

(27T) 2 



^,-K)n 33 (0,K^0) 

e 2 (k) 



A 2 (k) , E(k 
■ tanh 



1 



E 3 (k) 



2T 2TE 2 (k) 



cosh 



£7(k) 



2T 



(38) 



«a(v0) 2 } 



i- I «/'/ 



dn 

2^ 



d 2 K 



6(n, k)^^0(-o, -K)n^(o, k -> o) 



d 2 k 

(2n) 2 



- cosh 
2T 



2T 



(39) 



WFa(k)u_F/3(k) 



and the mixed ter m (p5[ ) does not contribute to the regular part within the used approximation, since II^O, K — > 
0) = 0. Evaluating (|38|) and ( p9| ) we performed the analytical continuation if2„ — » 57 + iO back to the real continuous 
frequencies, so that t is the real time. 

Using the nodal expansion ((2^) to calculate (38), (^9| ) and adding (|2^) we finally obtain the regular effective 
Lagrangian, C R such that /Jfikin = — * f dt f d 2 r£ n (i,r) for T< and ignoring the Landau terms: 

£* = —£dt0(t, r) + y (d t 6(t, r)) 2 - -(X78(t, r)) 2 , 
where the phase stiffness J = J^s and compressibility K = -Kd, s are 



(40) 



Jd = 



V<irv F v A 
24a 



In 2 bf„ 
T , 

27T 



1 



4a ^/ttvfva 



(41) 



The linear time derivative term in (^) is important for the description of vortex dynamics (see [^| and Rcfs. therein), 
but we omit it in what follows. We stress that the second, temperature dependent term in Jd follows from Eq. ( |39| ) 
which contains the derivative of the Fermi distribution, dnp(E)/dE. Its presence, as we mentioned in Sec. IV, makes 
the nodal approximation valid Since we consider only the low temperature T <C A^ region, we restrict ourselves 
by the values and va at T = 0, so that all temperature dependences appear to be linear. For higher temperatures 
it is necessary to take into account that Ad(T) and Va{T) are in fact decreasing functions of T. 

It is useful to compare the stiffness and compressibility from (41) with those parameters derived in fl5| for the 
continuum 2D model with s-wave pairing: 



Js = 



Am 



( 



\ 



\ 



dx- 



cosh 



^1 

AT 2 I 



K = — 

47T 



(42) 



First of all one can see that for low temperatures the superfluid stiffness in (^2|) does not contain any term which 
goes to zero more slowly than exp(— A s /T), while (f4l| ) has a term proportional to T. The origin of this difference is 
well-known and related to the presence of the nodal quasiparticles. Secondly, we can compare the values of the zero 
temperature superfluid stiffness which for the continuum translationally invariant system has to be equal to nj/4m, 
so that all carriers participate in the superfluid ground state ^5|. Since the presence of a lattice evidently breaks 
the continuum translational invariance the superfluid density at T — in the case of ( f4l| ) is less than nj /4m, as can 
be readily seen from ( p3[). Finally, since we consider here a neutral system it has the Berezinskii-Kosterlitz-Thouless 
(BKT) collective mode [{26|, D, 2 — v 2 K 2 with v = \J J/K. One can see that for the continuum s-wave case v = vp j\[2 
and 



'ttvfva 21n2at>i? 



(> 



(43) 
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for the d-wave model on lattice at T ~ 0. Eq. ([l3]) gives rather simple approximate expression for the BKT mode 
velocity for the lattice model of d-wave superconductor. Comparing in Fig. [2] the results obtained for v(T — 0) using 
Eq. ( [43| ) with the numerical computation without the nodal approximation, we can see that that even being very 
simple Eq. (0) predicts the correct behavior of v(T = 0). 

Following [Q] we estimate the upper values for the frequency Q, in (|3^) and the momentum K in ([39]). The 
"phase-only" effective action is appropriate for phase distortions whose energy is smaller than the condensation 
energy, E con d ~ iV(0)A 2 /2, where N(0) is the density of states. For the s-wave case ( ]42] ) this leads to the following 

restrictions: £1 < A s , vpK < A s and for the d-wave case (^lj): f2 < Ady/l/ao, v F K < Ad^fao, where an is the 
anisotropy of the Dirac spectrum defined in Introduction. 

VII. THE IMAGINARY PART OF THE LANDAU TERMS FOR 2D d-WAVE AND s-WAVE CASES 

The key values which are necessary for evaluation of the Landau terms are the differences E + — E- and n F (E + ) — 
uf(E-). Expanding in K, 

E(k + K/2)-E(k-K/2) = v g (k)K, (44) 

where the group velocity is given by 

v 9 (k) = V k £(k) = -L [£( k )v F + A(k)v A ] . (45) 

It is obvious that due to the gap k-dependence Eq. ( |44| ) differs from the s-wave case Q , where the difference is simply 

£(k + K/2)-^(k-K/2) = |^-v i rK (A(k)=A 5 ). (46) 

E{k) 



It is convenient to rewrite ( f44|) and ( |45| ) in terms of the nodal approximation described after Eq. ([25|) . In the vicinity 
of one of the nodes we have 

v s K = vpKiCOSLp + VAK 2 smp = P CQs(<p - ifr) , E(k ± K/2) < A d , (47) 

where the momentum K = (K\,K2) of 0-particle was also expressed in the nodal coordinate system ki, k2, so that 
P 1 = v F K x = Pcosip, P 2 = v A K 2 = PsimP and P = ^(P 1 ) 2 + (P 2 ) 2 = y/v'j,K'( +v 2 A K 2 . (We denoted the 
components of P as P 1 , P 2 to make them different from the node label Pj used in what follows.) 
The corresponding substitution in the integrals over K reads similarly to Eq. (B5ft 



(2tt) 2 ^ Jo 2ttv f va Jo 2i' y ' 



where P ma x is evidently related to the maximal value of K discussed at the end of Sec. VI , Finally, we can approximate 
the difference tif(E + ) — np(E-) as 

n F (E + ) - n F (E-) « ^v s K = ^Pco^ - *) (49) 
Having the differences ( pT| ) and (49) we can now derive the imaginary part for the Landau terms. In the subsequent 



subsections we consider all three terms of Eq. (|29j) for eZ-wave pairing and compare them with their s-wave counterparts. 
We would like to note that the Landau terms have also the real part which for the 3D s-wave case was consider in 
detail i n B . This real part consists of regular and irregular terms. The regular term was already taken into account 
in Sec. [Vj and the irregular term is not considered in this paper. 
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A. n^l{9Q 2 9} term 



Let us consider firstly the contribution from oS{0n^0} (see Eqs. 



continuation iQ n — ► Q + iO takes the form 
lm[i0Q^{m 2 e}] « - 



Oj) and (A2) for AJ), which after the analytical 



dttd 2 K , ,n 2 



-n,-K) 



d 2 k 1 
(2^)2 2 



4 

E 



dft 
2^ 



E 2 
2ttvfva 



Im- 



dn F {E) 



VgK + - 



"2V 



e(n,K) 



iO dE 1 
tt 2 



-fi,-K) 



pdp 




4 

E 



27TOf?;a 

dtt '■ Pr 

2^ 



d?2 ('El 

d^cos 2 fS(Pj cos(</3 — ipj) + fi) — 3 -fj cos(y 



d£ 



(50) 



27r^i7fA 



f2 




14 2 / 

-p2 COS 



2tt V ' ; 8 



sin 2 ipj 



-fi,-K) 



In 2 T 



7T WfWa 



e 1 



Here 0(x) is the step function and we used the integral 

pdp 1 



1 



In 2 1 



2nv F v A 2T cosll 2 V_ 
2T 



-T. 



TT VFVA 



(51) 



Integrating over tp we had to take into account that there are two points where (5-function contributes into the integral: 
cos<^ = —Q/Pj, simp = \ 1 — fl 2 /P 2 and cos<^ = —Q/Pj, simp = — ./l — Q 2 /P 2 . 



As we can see from the first equality in ( J5 0| ) the imaginary part develops when Q = E + — E- w v gK. This 
condition was interpreted in Q as the "Cerenkov" irradiation (absorption) condition for the process: "thermally 
excited ^-fluctuation quantum (phason) + quasiparticle < — ► phason + quasiparticle" (or phason being absorbed 
and scattering thermally excited quasiparticles). As was discussed in [[3), since definite energy is carried by the 
quasiparticles, this process is defined by their group velocity v 5 — dE(k)/dk as well as thermal and spin currents. 

It is in fact a coincidence that for the s-wave superconductor the directions of v g and vj? are the same, so that for 
the 2D s-wave superconductor using (46) instead of (Q) and ( ^5|) one can obtain 

d 2 k 2e 



lm[zf3Q^ n {0Q 2 0}} = - 



dQd 2 K T o 2 



f-n,-K) 



(2tt) 2 E 2 



Im 



— v F K + n + iO 

E 



ME) g 
dE E 



<md 2 K..^ m 2 



(— fi, — K) 



A s mc 
2tt 



dy 



dn 2\y\ 



'1-6 



i + y 

y 2 



dE^/T 



-.e 



(52) 



where following [0 we used the notations y = £/A. s and c = Q/vfK. Comparing ( |52| ) with the 3D case Q (in the 
notations of (J] the term we consider is related to the sum B^+Cl) one can notice that the only difference between the 
corresponding expressions is in the square root in the denominator of ( |52| ) . It obviously originates from the different 
measures for the angular integration in 2D and 3D where the extra sin ip is present. Comparing also ( |50|) and (|5^) one 
can see that their main analytical structure appears to be the same, viz. ~ Q 3 /Pj for the d-wave case and ~ Q /v F K 
for the s-wave case, so that the momentum variables Pj and v F K do not appear in the numerator. 

We note that the calculation of the ultrasonic attenuation a soun d(T, K) in d-wave superconductors results in the 
expression 



^sound 



d 2 k £ 2 (k) 



1 



(2tt) 2 E 2 (k) 



cosh' 



2T 



■S(y g K) . 



(53) 
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which has the same structure as (|50|) . This can be easily seen if one takes into account that for ultrasound frequency 
range f2 <C v F K, A^, so that the corresponding terms from Eq. ( |50|) can be simplified as follows 

S(E+ - E_ + tt)[n F (E-) - n F (E+)} = 5(E+ - E_ + tt)[n F {n + E+) - n F (E + )} ~ S15(E+ - E-) np }^ (54) 

dE 



leading to Eq. (|53|). This t he re ason why in the limit Cl/v F K — > the angular dependence of the Landau damping feCj ) 
which we consider in Sec. VIII would appear to be the same as the angular dependence of the ultrasonic attenuation 



n^' n {9K 2 e} term 



The contribution from nj^{0K 2 0} (see Eqs. (M\ and (|A2|) for A_) can be treated in the same manner 



d 2 k 1 
J2wf 2 
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E 

i=i 

Pn 
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E 



i 



pdp 

27TV F Va 

dn r Pm 

2^ Jo 

n 



2 P.J 



e 



i 

~ ~2 

e 2 + A 2 
e 2 

2ttv f va 



Im- 



v 9 K 



O + iO 



u_f q (k)w F/3 (k) dnF J^ v„K 



dE 



* 9( a,K,S^ 



0(-fi,-K) 



dipS(Pj cos(cp — ipj) 



2itv f va 

p 



2tt 



P 2 cos 2 ^ _ „ „,hi2 7 



(55) 



0(-O,-K) 



7T V F VA 



and for the s-wave case 



Im[0fi{^{0K 2 0}] 



dVLd 2 K 



dfld 2 K 
(2tt) 3 



0(fi, K) 



0, 



(2tt) e 



-0(fi,K) 



ft 2 



f-n,-K; 



A S TOC 

2tt 



-n,-K) 



dy 



d 2 fc 

■7 rrrlm^r- 

(2tt) 2 j_ 

E 



dn 2(1 + y 2 ) 3 / 2 



,2 dE 



2/ 2 M 



e 



-V Fa V F f3- 



vpK + + iO 



di(E) £ 
dE P 



(56) 

Let us now compare the expressions (^5|) and ( |56| ) . Evidently their analytical structure is different since for the s-wave 
case we have again that it is proportional to Cr/v F K, while for dwave pairing Pj enters the numerator, so that the 
main structure (excluding the measure of integration over Pj and the square root in the denominator) is ~ QPj . This 
difference originates from the angular integration in ( pq ) and ( |56| ) which is performed using the (5-function. Since the 
arguments of the (5-function for the d- and s-wave cases are different we have obtained two different answers. Indeed, 



for the s-wave case the argument is proportional to vj?K that coincides with the product (vfK 



2 dn F (E) 
' dE 



■v F K. outside 



the 5-function. (see Sec. after Eq. (|34|), where the origin of the term ~ (v^?K) 2 is discussed), so that the angular 
integration removes K from the numerator. The dwave case appears to be different since we have v s K inside the 



(5-function and (vj?K) 



, dn F {E) 
dE 



v g K outside, so that the angular integration leaves Pj in the numerator. Despite the 



fact that Q,Pj looks substantially less singular than fl/Pj, the first expression still remains nonanalytical near Pj ss 
because it is expressed in terms of VK 2 the coordinate representation of which is the non-local operator \J — V 2 . 

Thus physically the difference between the analytical form of the Landau damping in d and s-wave superconductors 
originates from the k-dependence of the gap A^(k) which in its turn makes the direction of the quasiparticle group 
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velocity v g different from the Fermi velocity vp. The last velocity, however, still enters the numerator of ( pq ) 
since, as was already mentioned in Sec. 0, it originates from the current-current correlator. The electrical current is 
proportional to vp, not v g because quasiparticles carry definite energy and spin, but do not carry definite charge. 
Therefore, the specific form of the Landau damping terms in a d-wave superconductor has the same physical origin as 
the source of the extra terms in the thermal and spin conductivities ^ which are related to the presence in v g both 
vf and va components. 

The s-wave Landau damping ([36|) itself can be again related to the 3D case j|, where it is expressed via the 
difference Bl — Cl- 



C. n£l{6Ktt6} term 
Finally we consider the contribution from fl < ^ n {9K.fl9} (sec Eqs. (^5|) and ([A3])): 

Im[»73fig>{0Kfi0}] « -\ J ^^6(Q,K)QK a e(-n,-K) 

d 2 k £ 1 ^ dn F (E) 

(2^ E lm v a K + n + i VFa (k) -^E~ V * K 

«E t£ f^p^ r^,K)^cos^,-K) 
j^J 2ir J 2nv F v A J 2ir V > 4 ^ V ' 

f p pdp f 2 * ^ dn F (E) n (57) 

x / / dip cos <pd(Fj cos((p — ipj) + II) — — — Pj cos(y> — ipj) y ' 

1 2ttvfv\ Jq dh 



7T V F VA 



2 , „ L 1^ 




_ cos ^eii- — 



and for the s-wave case 



M«mS{«n»H - -i / ^(n.Kjniwt-ii.-K) / 1 / f „^)| VfK 



^v f K + Q + iO 

E 



r*2 _ ' 



2 



9 

V 



Once again the difference between d- and s-wave cases is related to the different directions of v g and vp. Since this 
dtipiE*) 

time we had cos0vfK — ^— — v 5 K outside the (^-function, we could not get QPj as in Eq. (|55|) and obtained again 
dE 

~ Q 3 /Pj as in Eq. (p0[). Furthermore, because we have cos<£ instead of cos 2 tp as in Eq. (BfJ) we got only got only 
fl/ Pj cos ipj (compare with the square brackets in Eq. (|5Cj)) which after multiplying by IrP., gives the final result 
~ il 3 /Pj. It is interesting to note that the mixed term for d- (see Eq.(|57|)) and s-wave (see Eq.(j5|)) has the opposite 
with respect to other terms sign. This means that the mixed term describes the energy transfer from the quasiparticles 
to the phase excitations. It turns out, however, that the whole sum ( |59| ) of the Landau damping terms has the sign 
which corresponds to the phason damping. (We stress that the different sign in front of all final s-wave expressions 
and their d-wave counterparts are due to the explicit presence of the derivative of the Fermi distribution, which is 
negative, in the former expression.) 
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D. Final expression for the Landau damping term in the s-wave case 



Since the Landau damping terms for s-wave pairing all are ~ £1 /vfK, we can combine (52), ( pq ) and (|E 

1 f dVLd 2 K 



2tt 



dy dn 1 Q ( \y\ u \ (59) 



,i + y 2 dE 2 y2\y\y/TT7 \V^+y 2 



V y 2 

This expression coincides with the function H\ (c) introduced in pj except for the above mentioned difference between 
angular integration in 2D and 3D. Thus the simultaneous derivation of the s- and d-wave cases allowed us to be sure 
that all relevant terms were included and to see explicitly why the corresponding s- and d-wave terms behave so 
differently. 

E. Temperature and energy-momenta dependences of Landau Damping 

We now compare the conditions on the momenta and energies for the existence of the nonzero Landau terms in 
d- and s-wave cases and the temperature dependences of these terms. As one can see from Eqs. (j50|), ( j5q ) and ( |57| ) 
for a given node the imaginary parts develop when Pj > |f2| (or \Jv F K\ + v'^K$ > |0|). Thus in contrast to the 
s-wave case when the imaginary part is nonzero for all directions of the phason momentum K satisfying the condition 
— «f|K| < fl < v_f|K| (or |c| < 1, where c is defined after Eq. (p2|)) , for the d-wave it is sensitive to the direction of 
the phason momenta as shown in Fig. [|. As we will discuss in Sec. VIE for HTSC vf ^> Va, so that the directional 



anisotropy of the Landau damping becomes very strong. Indeed as can be seen from Fig. |^ for \CI\/vfK < 1 and 
an 1 the Landau damping would exist only in a narrow region of the momenta directions. Furthermore, if the 
projection K\ of the phason momentum K is not exactly zero (K \ k^), the cond ition |fi| < Pj can be well replaced 



by |fi| < We will discuss the validity of this approximation in Sec. VIIL It is important to stress that the 

sharp directional dependence discussed here is not related to the nodal approximation and follows only from the gap 
anisotropy and the fact that the difference of the Fermi distribution functions ( |49| ) which is present in Eqs. (^0|), ( |5q ) 
and ( prj ) has a very sharp k-dependence, so that in principle unbounded (for example, Uf(k) = k/m for the model 
with the quadratic dispersion law functions w s (k) and «F(k) in the corresponding integrals were replaced by their 
values at the Fermi surface. 

The final expressions ([50|), (|55|) and (|57| ) for the d-wave case are in fact simpler than their s-wave counterparts, J5^), 
( |56|) and ( |58|) (or the final expression (|59|)) because the latter expressions still have one integration which depends on 
the ratio |c| via the 0-function. Thus if |c| = \VL\/vfK — ► 1, only large y — £/A s contribute into the corresponding 
integrals. This circumstance and the opening of an isotropic gap A s make the contribution from the Landau terms 
very small comparing to the main terms. For d-wave pairing a relative contribution of the Landau terms does not 
depend on the ratio Q /Pj via p integration, so that all frequencies and momenta satisfying the conditions imposed 
by the 0-functions in ( |50| ) and ( |55| ) have the same temperature dependent weight. 

As in the case with the superfluid density the temperature dependence of the Landau terms is linear in T, but 
while the superfluid density is nonzero at T = 0, there is no Landau damping at T — 0. The linear T dependence and 
the absence of damping at T — are obviously related to the fact that for d-wave pairing there are only four gapless 
points on the Fermi surface. It is known, for example, that for a normal (nonsuperconducting) system the Lindhard 
function (or polarization bubble) 

tic, v\ f d2k -"*■(£-) / finA 

L ^ = J (27j H + -C-+» + *0 ' (60) 

has a nonzero imaginary part even at T = E3| if the Fermi surface remains ungapped since the derivative of the 
Fermi distribution at T — > becomes singular on the entire Fermi surface (or line in 2D). 



Finally, we note that there is also imaginary contribution even from the first, "superfluid" term in (A2) for A-. 
As we already mentioned, this term involves pair breaking which for s-wave pairing has the threshold energy, 2A S . 
For d-wave pairing this energy is zero for K = due to the presence of nodes. Nevertheless, if K ^ a finite energy 
Pj is necessary to create these excitations and the imaginary contribution is nonzero only for |f2| > Pj. Besides the 
analytical structure of this term is regular and it has a higher order than the terms considered here. Thus this term 
appears to be less important than the Landau terms we just considered. 
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VIII. THE APPROXIMATE FORM OF THE EFFECTIVE ACTION AND 0-PROPAGATOR 



Whereas the local nodal coordinate systems (Pj,tpj) are convenient to write down the corresponding contributions 
from each node, the final result should be presented in the global or laboratory coordinate system (K,<f>). It is 
convenient to measure the angle cf> from the vector k x , so that <j) — corresponds to the corner of the Fermi surface 
(see Fig. |l|) and the first node is at (j> — 7r/4. Thus the transformations from the global coordinate system into the 
local system related to the j-th node are 



cos Ipj 



F ' 



7T 71 , . 
-4 + 2C7-I 



2 • 2 
v A sin 



71 71 , , 

-4 + 2^'- 1 ))' 



sin 



v A K . 
■ sin 



(61) 



1,- 



,4. 



The estimates of show that in YBCO = vf/va — 14 and in BSCCO an — 19 (see Table |). Thus we can 
also use the inequality ap> 3> 1 in what follows. First of all this inequality implies that for j-th node we may assume 
that 



P 7 w Kv F 



(62) 



This approximation is, of course, valid only if K is not parallel to the Fermi surface (K !/( k 2 ). We note that 
this direction which is "dangerous" for j-th node is the nodal direction for the neighboring nodes. The size of the 
"dangerous" direction where Eq. ( |6^ ) becomes invalid can be estimated from the condition Kvf sin Acf> w Kva cos A(j> 
which gives A</> ~ a^ 1 . Since cud > 15 Eq. ( |62| ) is in fact w ell just ified outside the nodal regions. 

Using Eq. ( p2|) we can also rewrite the discussed in Sec. VII E inequality |f2| < Pj imposed by 8-functions which 
define the region where the Cerenkov condition can be satisfied as 



in I < v f k 



(63) 



One can easily see that the condition ( |63| ) is in fact equivalent to the condition |f2| < Vp we already mentioned. 
The effective action (H8h in the momentum representation in the global coordinate system can be written as 



with the regular (compare with Eqs. ( fl(i| ) and (|4l|)) 

K 2 ( y/irvFVA 2 In 2 up, 



F R (Q, K) 



4 V 6a 



-T 



71 VA 



1 

4 ay/nvpVA 



(64) 



(65) 



and Landau, .F L (fi, if, </>), parts. 

It is possible to write down F L (f2, if, <f>) substituting (|6l]) directly into ( p0| ) and j55|), but to have a more transparent 
expression we would like to consider a more simple case |fi| /vpK <C 1. This condition becomes equivalent to 
|f]|/Pj -C 1 due to the presence of 0-functions in (p0[), ( |55| ) and (57) which are cutting out the forbidden domains 
with |f2|/Pj > 1. Physically, the condition Q/vfK 1 is relevant if one, for instance, estimates the Landau damping 
for the BKT mode Q because it follows from (|f) that \Vt\/v F K = ^71/(60^ ) < 1. Thus we obtain 



. il 3 In 2 T 

' if 27T WA^F 
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(66) 
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(67) 
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\vpK ' 



cos 2 (0 - 7r/4) 



, ., e(|co8(0-7r/4)|--!^L 
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101 



sin 2 (0 - 7r/4) 
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(69) 



The functions fx, f% and f% are obtained from 
the assumption |fi|/Pj <C 1 replacing the square roots in 



from (50). Nevertheless, we kept the 0-functions which are present in 




respectively. Deriving (|67j) - (|69|) we used the 
1571) by 1. Furthermore, we kept only sin 2 ^- 
0|), (55) and (|57|) because they are essential 



in imposing the condition O/Pj. 

In Figs. |^ (a), || (a) and [j] (a) we show the functions fx, fi and f% which describe the intensity of Landau damping 
(/i and for the term ~ tt 3 /K and fo, for the term ~ VlK , respectively) as a function of the direction in the plane. 
Although the functions fx and fs describe the directional dependence for the same ~ £l 3 /K term, we do not combine 
these functions into the single function because the y o riginate from the different expressions. (We recall that fx 
originates from Eq. (50) and f$ originates from Eq. (p7|), respectively.) For the comparison in Figs. |] (b), [| (b) and 
P_[b) we show the directional dependences calculated by the direct substitution of Eq. ([n]) into Eqs. (p0"|), ( p5[ ) and 
(57) without making the approximations we just described. As one can see, the approximate representation (]6(|), 
(37) - (^) despite its relatively simple form (we have used only the terms ~ fl 3 /K and ~ flK) gives reasonably 
good expression for the Landau damping terms. The biggest discrepancies are seen between Fig. |] (a) and Fig. [| (b) 
because more approximatio ns w ere made to obtain fx term. 

As we mentioned in Sec. VII, the angular dependence described by fx (see Fig. ^] (a)) coincides with the angular 
dependence of the ultrasonic attenuation jl(J in the limit Q/vpK — > 0. Indeed the 0-functions from ( |67[ ) disappear 
when Q/vpK — > and Eq. ( |67| ) reduces to the corresponding equation from ]io| . 

We stress also that the analytical structure of the damping terms in ( |66| ) appears to be quite different from the 
dissipation introduced in ]27| into the "phase only" action basing on the low frequency conductivity measurements. 

Since fx + 2/ 3 and fi terms are odd functions of the frequency O, they integrate to zero in fltin- Nevertheless, the 
damping terms are manifest in the equation of motion and in the propagator of the BKT mode which is considered 
below. 

Comparing Figs. |^, || and ^, one can see that the functions fx and have more pronounced angular dependence 
than $2- Although as we explained above, the contribution from a given node is zero if for some K the condition ( |6^ ) 
is not satisfied, this condition can still simultaneously be satisfied for the neighboring nodes, so that the sum over all 
nodes in fx - fz is not necessarily zero. We note also that by the absolute value both fx (or fa) and f% terms are 
practically the same. This can be easily seen if we rewrite, for example, fx term as 



fi 3 In 2 T 
17 Tj 1 

K 7T V/\Vp 



n 



vpK 



vpK 



In 2 T ( n 



/: 

7T VA 



In 2 T 

aK— — fx 

b vp 



6m 



(70) 



where writing the last identity we explicitly used the ratio U,/vpK for the BKT mode at T = given above. 
In such a way following B we can write an approximate propagator of the BKT mode near T = as 



D?(fi,K) 



1 



4ay/irvFVA 



fl 2 - K 2 



/ TTVpVA 

\ 6 



2 In 2avp 



2ia-y(4>)TQK 



(71) 



with 



7 (0)= In 2 




(72) 



As discussed in [[| Eq. (fn]) has the form of a bosonic propagator with damping and its line-width has an explicit K 
dependence. In contrast to the s-wave case, the width depends not only on the absolute value of |K|, but also on the 
direction of K. The angular dependence for is shown in Fig. ^. Despite a simple form of the propagator (fnj), 
the corresponding effective wave operator (the transform of D7 1 to space and time) has a nonlocal form in coordinate 
space due to the presence of the damping term which depends H on |K|. 
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Comparing the second term in parentheses of ( |7lD with the damping term one can see that they have the same 
order of magnitude showing that even for the low temperature region the Landau damping becomes important and 
its magnitude is comparable with term which describes the linear low temperature decrease in the superfluid stiffness. 

We stress also the difference between the Landau damping representation in Eqs. ( |7l|) and (66). While the propa- 
gator^^) relies on the particular dispersion law for the BKT mode with the approximate v given by Eq. (]43|), our 
Eq. (|66j) along with Eqs. (|67| ) - ((6^) present a rather general representation for the Landau damping terms derivation 
of which did not rely on any particular dispersion law for the phase excitations. Thus it can be, in principle, used to 
describe the damping of the plasmons in charged superconductor. The only assumption we made is that n/vpK <C 1, 
which was used to derive more simple an d tra nsparent representation for the Landau damping. For the more general 
case one should use the results from Sec. VIL 



IX. CONCLUDING REMARKS 

It is very important to note that the collective phase excitations described by the propagator similar to ( f7l| ) can and 
have indeed been studied experimentally. Indeed the measurements of the order parameter dynamical structure factor 
in the dirty Al films allowed to extract the dispersion relation of the corresponding Carlson-Goldman mode and to 
investigate its temperature dependence |28| . It is important to stress that since the real systems are charged this mode 
appears to be different from the sound-like Bogolyubov- Anderson mode which exists in a neutral superfluid Fermi 
liquid. However, as discussed in [^2| (see also ]§| and Refs. therein) the discovery of Carlson-Goldman collective mode 
p8j overcame the widespread opinion that the Bogolyubov-Anderson collective oscillations predicted for a neutral 
superfluid should be forced to the frequencies on the order of the plasma frequency. It appears that if under certain 
conditions the frequency of the phase excitations is sufficiently low, the normal electron fluid may screen completely 
the associated electric field, thereby making the theory of uncharged superconductor applicable [M . The conditions 
for screening and occurrence of the Carlson-Goldman mode in s-wave superconductors are favorable only if the systems 
are dirty (to suppress the Landau damping) and T <T C . As recently claimed in Q in d-wave superconductors the 
conditions appear to much less strict, so that the Carlson-Goldman mode may be observed in the clean systems for 
T down to 0.2T C . Now there are some new questions which would be interesting to address from both experimental 
and theoretical points of view. 

First of all, there is a general question whether the Carlson-Goldman mode which is the equivalent of the Goldstone 
mode for the charged system can be observed experimentally. There are also a lot of theoretical questions which study 
would make this kind of experiment possible. Leaving aside the specific of the Carlson-Goldman mode studied in Q 
we mention some of them which are directly related to the present work. 

1. The Carlson-Goldman experiment [^8| measures so called dynamical structure factor. This factor has a peak 
associated with the phase excitations the width of which is primarily controlled by the Landau damping. Therefore 
our result that the intensity of Landau damping has a strong directional dependence should be taken into account 
in the calculation of the structure factor. In particular, the lowest for the observation of the Carlson-Goldman mode 
temperature 0.2T C in J^] is obtained for the nodal direction. As can be seen from, for example, FigJJ and ||, the 
corresponding Landau damping terms become maximal in the vicinity of this direction which should likely result in 
widening of the corresponding peak. If this widening is too big that the peak may become unobservable, this would 
again suggest that the dirty samples should be used for the experiment. 

2. Taking into account the lattice effects and the dominance of the nodal excitations we have obtained that the velocity 
of the phase excitations at T — > is given by ui at = «f V /&Q-D (see Eq. ( f43| ) and the discussion after Eq. (p3|)). 
This expression appears to be quite different from the well-known continuum result, w CO nt = Vf/v2. The estimates 
of the velocities w C ont and i>i a t obtained from the values vp and ajj [p"l| are presented in Table |[ It is seen from these 
estimates that the difference between the description proceeding from continuum and lattice models can be quite 
different. However the estimates of the minimal for the observation of the Carlson-Goldman mode temperature in || 
are based on the continuum expression w con t = vp/y/2. Thus it would be interesting to reconsider these estimates for 
the lattice case. 

3. Finally, as pointed out in II if the plasmon is at finite frequency at K — ► 0, the Landau damping does not occur 
since when Q is finite and K — > it is impossible to satisfy Cerenkov condition discussed above. However, a very small 
value of the plasma frequency in HTSC suggests that the Landau damping may still be relevant when K is nonzero 
and the Cerenkov condition can be satisfied. Thus the corresponding damping terms should be included in the "phase 
only" actions which are used to describe plasma excitations in <i-wave superconductors. Predicted anisotropy of the 
Landau damping would result in the damping anisotropy for the plasma excitations. 

To summarize, we have considered the "phase only" effective action for so called neutral (or uncharged) fermionic 
superfluid with d- and s-wave pairing in 2D. When the damping terms are included into this action, it turns out 
non-local in coordinate space and its analytical structure for the d-wave case is very different from the s-wave case. 
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To consider a charged superfluid it is necessary to combine the approaches used in the present paper and in |g, |9J to 
consider the Landau damping in the presence of Coulomb interaction. 
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APPENDIX A 



Substituting @ into (f 
diagonal terms we obtain 



and (B4), using the definitions for m-j and IL^- and evaluating the matrix traces for the 



T E 



l= — oo 



d 2 k 2[{iuJi + iVL n )iuji + £(k - K/2)£(k + K/2) ± A(k - K/2)A(k + K/2)]V±(k) 



(Al) 



(2tt) 2 [(iui + in n ) 2 - £ 2 (k + K/2) - A 2 (k + K/2)][(^/) 2 - £ 2 (k - K/2) - A 2 (k - K/2)] ' 



V±(k) 



v Fa (k)v FI 3(k), 
1, 



Performing in ( |Al[ ) the summation over Matsubara frequencies and simplifying the result we arrive at 
A± = 



d z k ri 



E-E 4 



1 



1 



1 ( e-e + ±A_A H 

2 V E-E+ 



1 



E++E-+ ifl n E+ + E-- ifi 
1 



[l-n F (E_)-n F (E + )] 



(A2) 



E+-E-+ ifl n E+-E-- iVL r , 



n F (E_)-n F (E + )]}V±(k) 



where £± = f (k ± K/2), E± = E(k ± K/2) and A± = A(k ± K/2). The second term in Eq. (|A2j) can be further 
simplified if one notices that the term with 1/{E + — E- — ifl n ) transforms into the term (—1)/(E + — £L + iCl n ) when 
k — > k. Note that the A± terms do coincide with the corresponding terms in [Q. 
For the mixed term p5l) one obtains 
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-k and took into account th at U ffk) also changes sign under 
IlQ 3 (ii7„,K). Eq. (A3) is also in agreement with the 



To get the last identity we again replaced k - 
this transformation. One can also check that Il3 (if2 n ,K) 
corresponding term in Q denoted as D. 

The first and second terms in (A2) and (A3) have a clear physical interpretation (24]. The first term gives the 
contribution from "superfluid" electrons. The second term gives the contribution of the thermally excited quasipar- 
ticles (i.e. "normal" fluid component). The essential physical difference between the two terms is that the superfluid 
term involves creation of two quasiparticles, with the minimum excitation energy in the s-wave case being 2A S . (For 
the (i-wave case the minimum energy turns out to be zero in the four nodal points, see the discussion at the end of 

Sec. jyg .) 

On the other hand, the normal fluid term involves scattering of the quasiparticles already present and the excitation 
energy in this case can be arbitrary small, as in the normal metal, independently whether we are in the vicinity of the 
node or not. Of course, the number of the quasiparticles participating in this scattering depends on the value of the 
gap and drastically increases if the gap is equal to zero in some p oin ts. As we shall see, the Landau damping terms 
originate from the second, normal fluid term. Note that our Eqs. ( A2) and (|X|) (as well as Eqs. (|30|), (34) and (|3~F 
are suitable for studying both s- and d-wave cases. 
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TABLE I: The Fermi velocities Vf and anisotropies of the Dirac spectrum, an of YBa2Cu3C>6.9 and Bi2Sr2CaCu208 at optimal 
doping from 111). The velocities of the phase excitations for the continuum, v corLt and lattice, «i at limits at T = 0. 





Vf, x 10 7 cm/s 




Ucont, xlO 7 sm/s 


fiat, xlO 6 sm/c 


YBCO 


2.5 


14 


1.77 


4.8 


BSCCO 


2.5 


19 


1.77 


4.2 




FIG. 1: For fi = the Fermi surface (dotted line) represents the points where s(k) = — 2t(cos k x +cos k y ) = (in the units where 
the lattice constant a — 1). There are four nodes centered at (±7r/2, ±7r/2) around which the energy spectrum is linearized. 
The corresponding nodal sub-zones (see Eq. (pB[)) are called BZ j with j = 1, . . . ,4. For qd C 1 the Landau damping develops 
only if the direction of the phason momentum K is within one of the "cones" . The size of these cones is dependent on the ratio 
c = Q/vfK, so that the cones shrink as |c| — > 1 (see Eqs. ([^) - (|69|)). 
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FIG. 2: The dependence of the Goldstone mode velocity, v(T = 0) on the amplitude of the gap, Ad- Solid line is the result of 
numerical calculation with £(k) = — 2(cos k x + cos k y ), A(k) = Ad/2(cos k x — cos k y ) (we put t = a = 1, so that A^ is expressed 
in units of t). Thin line is obtained using Eq. (^3|). The anisotropy of the Dirac spectrum for this case is qd = 4/Aa, so that 
A d = 0.2 corresponds to qd = 20. 




a given only excitations with the momenta K which are outside the ellipse (in the shaded region) can contribute into the 
Landau damping. For a fixed ratio |Q|/i;f|K| the presence of the Landau damping depends not only on |K|, but also on its 
direction. For example, for a vector having the length of the vector K shown in the figure the Landau damping is possible only 
if its direction is sufficiently close to the direction of ki normal to the Fermi surface. 



21 



2 
1.5 
fi 1 
0.5 



0.5 



*7l" 



1.5 



2 
1.5 
fi 1 
0.5 






* nf 



FIG. 4: (a) The angular dependence, /i y-^j^,<pj given by Eq. (|^) of the Landau term ~ Q a /K for q d = 20 and \Q\/vpK 
0.1. (b) The same function calculated directly from Eq. (|5o|). 




FIG. 5: (a) The angular dependence, fa y V ^ K , 4>J given by Eq. ( |6q ) of the Landau term ~ Q.K. The parameters q_d and 
Ifil/i^F-K" are the same as in Fig. H. (b) The same function calculated directly from Eq. (pq). 
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